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1 Introduction 



Lattice Bose-gas models were invented as an alternative way to understand continuous inter- 
acting boson systems including liquid Helium, see |MMj and a very complete review [H]- But 
recent experiments with cold bosons in traps of three-dimensional optical lattice potentials show 
that lattice models are also relevant for describing the experimentally observed Mott insulator- 
superfluid (or condensate) phase transition |G-Bj . In |BDj and then in |A-Yj , this phenomenon 
was analyzed rigorously in the framework of the so-called Bose-Huhhard model. 

The aim of the present paper is to study a disordered Bose- Hubbard model and in particular 
the influence of the single-site potential randomness on the Bose- Einstein condensate (BEG). 

Notice that the first attempts to understand this influence go back to |KLlj . |KL2j and |LSj 
for continuous Perfect Bose-Gases (PBG) in a random potential of impurities. For the rigorous 
solution of this problem see |L-Zj . One of the principal result of |L-Zj is that the randomness 
enhances the BEG. For example, the one- dimensional PBG has no BEG because of the high 
value of the one-particle density of states in the vicinity of the bottom of the spectrum above 
the ground state, making the integral for the critical particle density infinite. The presence of a 
non-negative homogeneous ergodic random potential modifies the one-particle density of states 
(due to the Lifshitz tail) in such a way that the integral for the critical density becomes finite. 
Hence, the one- dimensional PBG with random potential does manifest BEG. The nature of this 
BEG is close to what is known as the " Bose-glass" since it may be localized by the random 
potential |LZj . This is of interest for experiments with liquid ^ife in random environments like 
Aerogel and Vycor glass, |F-Fj . |KTj . 

On the other hand, the nature and behaviour of the lattice BEG may be quite different. 
First of all, the lattice Laplacian and the Bose-Hubbard interaction produce a coexistence of 
the BEG [superfluidity) and the Mott insulating phase as well as domains of incompressibility, 
see e.g. |F-Fj . |K-Gj . Adding disorder makes the corresponding models much more compli- 
cated. The physical arguments |F-Fj . |K-Gj show that the randomness may suppress the BEG 
(superfluidity) as well as the Mott phase in favour of the localized Bose-glass phase, but this 
is very sensitive to the choice of the random distribution. 

Since there are very few rigorous results about the BEG in disordered systems, we consider 
here a single-site random version of the lattice Infinite-Range Hopping (IRH) Bose-Hubbard 
model, which in non-random case has recently been studied in detail for all temperatures and 
chemical potentials in |BDj . 

This paper is organized as follows. In Section 2, we define the lattice Laplacian for finite- 
and infinite-range hopping and recall the results about BEG for the free lattice Bose-gas. We 
then introduce random single-site and on-set particle interaction potentials and state our main 
result about the existence of and an explicit formula for the pressure for the IRH Bose-Hubbard 
model with these type of randomness. We outline the proof of the main theorem using the 
approximating Hamiltonian method. 

In Section 3 we consider the pressure for extremal cases of hard-core and perfect bosons. 
We show that they are the limits of the IRH Bose-Hubbard model pressure when the on-site 
particle interaction tends respectively to -|-oo and to 0. 

In Section 4, we analyse the phase diagram in the case of a non-random on-site particle 
interaction and random single-site external potential. We distinguish a number of different 
cases. We start with perfect bosons and show that the randomness enhances BEG in this 
case, see Sect.4.1. This is no longer true for interacting bosons. We study in Sect. 4. 2 the phase 
diagram first for Bernoulli single-site potential and then for trinomial and multinomial discrete 
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distributions. 

In the case of a Bernoulli distribution and hard-core bosons (infinite on-set repulsion) we 
showe that in addition to the complete BEC suppression at extremal allowed densities p = 
and p = 1 there is a new point p = 1 —p, where p = Pt {potential ^ 0}. We prove that for finite 
on-site repulsion the suppression of BEC at integer^ and also for fractional values of densities 
p = n — p ,n = 1,2,... persists, if the Bernoulli potential amplitude is large enough. In fact 
we find that increasing the Bernoulli potential amplitude (disorder) decreases the critical BEC 
temperature in the vicinity of fractional values of densities but increases it for integer values of 
density. A similar phenomenon occurs also for equiprobable trinomial distributions, but now for 
densities p = n/3. Our numerical calculations demonstrate that it should be true for a general 
multinomial distribution. 

For illustration of a continuous distribution we study a homogenous distribution with com- 
pact support. Then for hard-core bosons we prove that the complete BEC suppression occurs 
only at extremal allowed densities p = and p = 1, with the trace of suppressions only at 
integer values of densities for a finite on-site repulsion. In particular we show that the critical 
BEC temperature gets lower, when one switches on disorder for (a small) on-site interaction, 
whereas it gets higher for perfect bosons. For large values of on-site interaction the picture is 
similar to the discrete distributions: increasing of disorder increases the critical BEC temper- 
ature in the vicinity of integer values of density but increases it for the complimentary values 
of density. 

In Section 5 we summarize and discuss our results. 

2 Model and Main Theorem 

For simplicity we shall consider the Bose-Hubbard model only with periodic boundary condi- 
tions. So let A := {x E Z'^ : —La/2 < < L^/l, a = 1, . . . , c?} be a bounded rectangular 
domain of the cubic lattice TJ^ wrapped onto a torus. Then the set A* := {q^ = 27m/ : 
n = 0,±1,±2,... ± {La/2 — l),La/2, a = l,2,...d} is dual to A with respect to Fourier 
transformation on the domain A = Li x L2 x . . . x of volume |A| = V. 

The standard one-particle Hilbert space for the set A can be taken as f)(A) := with the 
canonical basis {e^.}^.gA, i.e. e^iy) = 5x,y. Then for any element u = ^^.^a ^^-^^^ ^ ^i-^) 
one-particle kinetic- energy (hopping) operator is defined by 

(tAM)(x) ■.= J2ttyiuix)-uiy)) = ^t^,j^(M, -«,), (2.1) 

where 

= 7 E ^"^^^"^^"'^ ' (2.2) 
(?ga* 

is the periodic extension in domain A of a symmetric, translation invariant and positive- definite 
matrix, i.e. 

4 = E</'"^0- (2.3) 

2/6A 

Notice that functions {{eq){y) := e^''^ / \/V}q^A* also form a basis in f)(A), i.e. for any u E f)(A) 

one has u = J2geA* '^1^1 ■ 
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Let '■= 5^s(l)(A)) be the boson Fock space over f)(A). For any / G [)(A)) we can associate 
in this space the creation and annihilation operators 

a*{f) ■.= J2a*{y)f{y) , a{f) := $^a(2/)r(y) . (2.4) 

3/eA j/SA 

Let a*, and a*, a, be the boson creation and annihilation operators corresponding respec- 
tively to the basis elements Cx and Cg, satisfying the lattice Canonical Commutation Relations: 
[ax,a*] = Sx,y and [ag,a*] = Sg^p. Then = 0*0^ is the one-site number operator, and 

NA:='^nx= '^a*ag, (2.5) 

is the total number operator. 

The second quantization of the hopping operator ()2.1|) in gives the free boson Hamilton 
of the form 

Ta := ^<(tAa)x = \Y1 ^^y^^*^ ~ "P*^"^ ~ ^y") = ^ito-iq)a*gdg. (2.6) 

X 6A x,y£K qSA* 

If hopping is allowed only between the nearest neighbor (n.n.) sites with equal probabilities, 
then = —A corresponds to minus the lattice Laplacian, i.e. 

d 

^xy = XI (^^+1- .2/ + '^^-1^ .2/)' (2-7) 

where (a; ± Iq,)/? = X/j ± (5q,/3. In this case the one-particle hopping operator spectrum is 

d 

e(g) := (to - ig) = 4 sin^(g„/2) > , g G A* , (2.8) 



a=l 



with eigenf unctions {egjggA*- 

It is known that the lattice free Bose-gas ()2.6p with n.n. hopping manifests the zero-mode 
BEC when d > 2, since the spectral density of states Md^de) corresponding to ()2.7p is small 
enough to make the critical particle density pl^'^^{P) bounded for a given temperature : 

P^nnl/^) -limlimi V , ,\ = ttt^ / A -^r^^ — 7 (2.9) 



Ji8+ e'" — i 



Here limA stands for the thermodynamic limit A f Z'^, by B"^ := [— 7r,7r]'^ we denote the first 
Brillouin zone and the density of states N'd{de) = {cd£^'^^'^~^^ + o{e^'^^'^~^^)}de for small e. 
A similar result is true for the infinite-range (i.r.) hopping Laplacian: 

tty = l^a-S.,y) , x,yeA. (2.10) 
By ()2.10|) the one-particle spectrum in this case takes the form: 

e(g) := (to - ig) = {1 - 5g,o) > , q E A* . (2.11) 
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Therefore, it has a gap: 

hme(g) = 1 ^ e(0) = , (2.12) 

and allowed values of the chemical potential are still n < 0. Since the density of states is simply 
zero in the gap, and |A*| = V^IS"^], we have Md{de) = 6{e — l)de. Therefore, the critical particle 
density has a bounded value: 

'■t-.l/S) ~/^/^, = ^ < oo , (2.3) 

for any dimensions. The latter implies a zero- mode BEC for densities p > p{^^ij. {(3). 

The problem of existence of BEC gets much less obvious if one takes into account the boson 
interaction. This is even the case for the simplest on-site repulsive interaction 

i/A:=TA + A ^n,(n, -1) , A>0 , (2.14) 

known as the Bose-Hubbard model. (Notice that attraction: A < makes this model unstable, 
see jn] for discussion of other cases.) 

Remark 2.1 Concerning the model ^2.14\ ) ^he best rigorous results so far are: 

- a proof of BEC for the n.n. lattice Laplacian and the hard-core boson repulsion: A = +oo, by 
|K-Sj for the case of the half-filled lattice, see also |ABj .• 

- a recent exact solution of the IRH Bose-Hubbard model \2.1(^) . \2.14^ for any X > by |BDj . 

The aim of the the present paper is to study a disordered IRH Bose-Hubbard model. Let 
[Q, S,P) be a probability space. We define our basic model by the random Hamiltonian: 

= ^ E « - - + E^>-('^- - 1) + (2.15) 

x,j/6A xeA xeA 

where parameters {X'^ > 0}^^^d and {e'^ G M}}^^^d, for u E fl, are real-valued random fields 
on Z'', which we suppose to be stationary and ergodic. We denote by 

p%{(3,p) :=p[Ht]i/3,p) := ±^Tr^^exp{-/3{HX-pNA)} (2.16) 

the grand canonical pressure of the system ()2.15j) for given temperature (3~^ and chemical 
potential /i. For non-random parameters Al^ = A > and e'^ = e = the model ()2.15|) was 
considered in |BDj . 

Our main theorem is a formula for the pressure of this model given some general regularity 
conditions on the random parameters involved in the Hamiltonian (j2.15p . 

Theorem 2.1 Let the stationary, ergodic random fields {AlJ^j^g^d and {st}^^id be such that: 

Amin := inf A^ > , emin := inf > -oo. (2.17) 

Then for almost all u G VL, i.e., almost sure (a.s.), there exists a non-random thermodynamic 
limit of the pressure \2.1(^) : 

a.s.-\lmpl{(3,p)=p{P,^^), (2.18) 

A 



such that 

p{f3,fi)= (2.19) 
sup {-r^ + P~^E {lnTr(5^)^ exp/? [(// - - l)n^ - A;;;n^.(n^. - 1) + r(a* + ctx)]}} , 

r>0 

where E (■) zs expectation with respect to the measure P. 
Proof: Let 

^OA := E - 1) + E(^^ + 1)^- • (2.20) 
Then by definitions ()2.4|) the Hamiltonian (j2.15p takes the form 

HX = Ta + J2 A>xK - 1) + E = -^S«o + i^o^A • (2.21) 

xeA xeA 

Since conditions ()2.17|) imply the estimate from below: 

HX > -a*ao + Na + A^i„ E '^-('^- " 1) + ^^iniVA (2.22) 

xeA 

> ^iVl + - A^i„)iVA , 

the Hamiltonian ()2.21|) is superstable. Thus, the pressure in fl2.18|) is defined for all /i G M^. 
Following jB-Tj . we introduce a similar Hamiltonian with sources: 

HX{u) := HX - VviVaQ + z/a*) , z/ G C , (2.23) 

and the corresponding approximating Hamiltonian: 

HX{z, v) := H^j,{z) - Vv{Vao + ua*,) , (2.24) 

where 

H^^iz) := H^j^ + V\zf -Vv{zao + zal) , z e C. (2.25) 

Then 

HXii^) - HXiz, u) = -(ao - zy/vr{a, - zy/V), (2.26) 
and by virtue of the Bogoliubov convexity inequality one gets the estimates: 

< p [HX{u)] - p [HXiz, u)] <^({ao- zWyiao - zW)) (2.27) 

for each realization G f2. Here {—)ii^(y\^ '■= (/^' /^) denotes the grand-canonical 

quantum Gibbs state with Hamiltonian ()2.23p . and from now on we systematically omit the 
arguments {f3,fi)- If we choose in the right-hand side of ()2.27p 

\ («o)hxM ' (2-28) 



then ()2.27|) implies the following estimate for each u & Q: 



< p [HXiu)] - snpp[HX{z, u)] < 1 {6al 5a,)^^,^ , (2.29) 



where we denote 

Sao := ao - {ao) jj^^^y (2.30) 

Since (|2.5|) implies the estimates: 

-Wiyao + ual) >-\iyf a^ao -V > -{u^Na - V, (2.31) 
by virtue of ()2.22|) and ()2.3H) the Hamiltonian with sources ()2.23p is also superstahle: 

Ht{iy) > + (£„,i„ - A„,i„ - \uf) Na-V , (2.32) 

uniformly in u; G and in \u\ < Cq, for a fixed Cq > 0. The superstability ()2.32j) implies that 
there is a monotonous nondecreasing function M := M{j3,fi) > of yU G M^, such that for any 
u; G n we have the bounds: 
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< {Na/V)^.^,^ (A/i) = d,p[HXiu)] (/3,/i) < M^ (2.33) 

and 

\zA,^if3,fi;iy)f <M^ (2.34) 
for the maximizer 2A,a;(i^) := ZA^i^{P,fi; u) in ()2.29p : 

P [^a(^a,c.(/5, /U; z^), J^)] iP, /i) := supp Kiz, z/)] (/3, /i) , (2.35) 

zee 

uniform in |z/| < Cq. Notice that the maximizer satisfies the equation: 

ZA,M = duP [H^izA^M, ^)] = (ao/W) . (2.36) 

Moreover, by the same line of reasoning as in |ZBj . Ch.4 (see also jBDj ) one gets that for 
|z^| < Co there are some u = u{M) > and w = w{M) > such that 

{6al 6ao) < {u + w(5a* , 6ao)HX{u) } , (2.37) 

where 

{5al , Sao)Hxi,) = p-'d, d^P [H^a{^)] . (2.38) 
Then the estimates ^TT^ and ^TTT) imply: 

< p [Hl{u)] - p [Hl{zA,M. ^)] < ^ + , ^«o)/fXH} • (2-39) 

Following |PSj we define in the Hilbert space L^({(Rei/, Imz/) : \v\ < Cq}) the Dirichlet 
self-adjoint extension Ly of the operator 

Lv := I -w{f3Vy^d,d^ . (2.40) 

Here 49j, djj = A coincides with the two-dimensional Laplacian operator in variables (Rei^, Imz/). 
The operator Ly is invertible and Ly^ has the kernel (^-^y^j {Green function), and 



^I/y^j {iy,i^') = for \iy\ = Co, or |z/'| = Cq, by the Dirichlet boundary condition. Since the 

semigroup |exp [—t{Lv — is positivity preserving, the same property is true for the 

operator Ly^, see e.g. |RS2j . Ch.X.4. 

Now, let p^u) := p [H'^{h')] and Po{i') := p [H^{z\^^{h'), u)]. Since Ly^ is positivity preserv- 
ing, then (I23ni)- (12301) imply 

(V(Po + V^))M>pW, (2.41) 
and by consequence the estimates 

< p [HX{u)] - p m{zA,M, 1^)] < [Ly\p, + u/V)) (z/) - p,{u) 

^Yu^u'){p,{u')-p,{p)} + u/V, (2.42) 



< / du' (Ly 

\y'\<Co 



where we used that du' {^v^ ~ 1^1 ^ ^o- By virtue of ()2.34p and ()2.36p we 
obtain for the integral in the right-hand side of ()4.4()|1 the estimate: 

/ du' fiy') (z/, u') {po{u') - po{u)} < 2M [ du' fiy') {u, u') \u' - u\ = ly. (2.43) 

J\i''\<Co ^ ^ J\v'\<Co ^ ^ 

After change of variables to = u\/V , we get 

= ^ / di' (L-yU) (e,0 W - < ^. (2.44) 
Here we used that in the Green function is known explicitly: 



where the Bessel function Ko{x) ~ A/'7r/2xexp(— x) decays exponentially fast for large x > 0. 
Therefore, and (jCTj) imply 

< p [HXiu]] - p [HXizA,^^), u)] < Oil/V) , (2.46) 

for all u E fl, any P > 0, ^ E M} and \u\ < Cq. 

Notice that by definitions ()2.20|) and ()2.25p for any 2;, z/ G C we get: 

pl^^pAP, /x; z, u) := p [HXiz, u)] iP, /i) = - \zf + (2.47) 
+ ^ 1^ Trj, ^ 1)'"- - ^>-(^- - 1) + (^ + + (^ + • 

Then ergodicity of the random fields {A|^}^g^d and {s'^}x^id implies the existence of the a.s. 
limit: 

P appr(/5, /i; 2, Z^) = a-S. - limp5(_ appr(/^' /^^ 2, z/) = - \z\^ + (2.48) 

+/3-iE{lnTr^^ exp/? [(/i - - l)n^ - AX(n^ - 1) + (2 + u)al + (z + z7)a^]} , 
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i.e., the self- averaging |PFj of the hmiting approximating pressure p'^ppj.(/?, /x; ^, i^). 
Now we put the source u —>■ and we make the canonical (gauge) transformation: 

■= a^e'^'^'. (2.49) 

Since Hamiltonian (j2.25j) is invariant with respect of this transformation, we get that z = \z\ := 
r and fcf. dTTTfl ): 

PA,appr(/5, /i; r) := PI,,JP, ^;z = r,u = 0)=p [H^ir, 0)] ^^) = -r' + (2.50) 



^ ^ In Tr^^ exp p [(/i - - l)n^. - A>^.(n^ - 1) + r(5; + 



Therefore, without source the maximizers in (|2.35|) can be defined only up to a phase and their 
moduli satisfy the equation: 



^ = ^ E + ~<^HXir,o) =■■ air), (2.51) 

where 



2^ . 



(r) '■— {Clx + '^x) Hf^{rfl) ~ 

^ Trg, {{a^ + 5*) exp/j [(/x - - l)n.^ - X%n^{n^ - 1) + r(a: + ~a^)]} 
Trj. exp /3 [(/x - - l)n^ - A^n^(n^ - 1) + r(a* + S^)] 



(2.52) 



When r = 0, the approximating Hamiltonian ()2.25|) is invariant with respect to canonical gauge 
transformations U^^cixU* = a^je*'^ for any (p. This implies ^'^{r = 0) = 0. Hence, equation ()2.5ip 
always has a trivial solution r = and , moreover, by ()2.34|1 any nontrivial solution < M. 
Finally, differentiating ()2.52|) with respect to r we obtain: 

< drCir) < R, (2.53) 

where, by the superstability (j2.32j) . the upper bound R is finite uniformly in uj,r,x. Hence, 
-2M < drP%^pp^{(3,fi]r) < 2RM for r e [0,M]. By consequence the limit (ITiHD implies the 
uniform a.s. convergence of the sequence {PA,appr(/^' '^)}a ^ ^ [0,^]: 

Pappr(/?, /i; r) = a.s. - limp;;(_^pp,(/3, /i; r) = (2.54) 
= -r^ + /?-iE{lnTr5, exp/5 [(// - - l)n,. - A>,(n, - 1) + r{al + a,)]} , 

Therefore, 

a.s. - limsupp^ (/3,/i;r) = suppappr(/?, /i; r). (2.55) 

A r>0 r>0 

Together with and (jTiH|l . the hmit (E3K|l proves the assertions and (|TTn|l of the 

theorem. □ 

Remark 2.2 The function C,x{f) is increasing in r by virtue of \2.5'J) . Moreover, it has also 
been suggested that for any x G Z'^ and uj G VL, the function r i— > ^^(r) is concave, see jBDj 
for discussion of this conjecture. This implies that the nontrivial solution of equation i2.51}} 
is unique. Notice that homogeneity and ergodicity of the random field random field {£x}x£Z'i 
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implies the same for the random field {^x}xez'' defined by i2.5^) . Therefore, equation \2.51\) in 
the thermodynamic limit takes the form: 

r = a.s. - limilir) = ^ E (C=o(^)) =■ fir), (2.56) 

expressing a self- averaging property of the order parameter r , see [EE]- Since the expectation in 
i2.5t^} preserves convexity, solution of the limit equation \2.5b]) should he also unique. There- 
fore, the sequence of maximizers {r%}^ with P = 1 has a unique accumulation point in the 
interval [0,M]. Moreover, if r'^ is the unique solution of equation i2.51]) . then 

a.s.-\imr% = r{p,fj.), (2.57) 

A 

where r{j3,fi) denotes the unique solution of equation \2.5f^) . 

Proof: Since Amin > 0, by superstability we get < M, see fj2.34|) . i.e. 

< lim inf < lim sup < M, (2.58) 

for any u E ^. Now suppose that there exists f2> with P(f2>) > and a subsequence 
{rln]n>i^^ ^ ^> ^^^^^ 

hmrX, =<>r(/3,/i), a; G ^]>. (2.59) 

n.— >oo 

Then, by virtue of (jSIHl), (I23SI), (f^3n|) and we get: 

- R |rt - < I < rt^ = a,Sr: + rt - C) < ei(rr) + R \rl^ - | . (2.60) 

These estimates, together with the hmit ()2.59|) and a.s. -convergence of ^a„('") /(^) ^ 
imply 

< = /(0>^(A/i), (2.61) 

for any uj G f2> with P(f2>) > 0, which is impossible by uniqueness of solution of ()2.56|) . 
Similarly one excludes the hypothesis rif < r{j3,fi), which proves ()2.57|) . □ 



3 Limiting Hamiltonians 
3.1 Limit of Hard-Core Bosons 

The hard-core (h.c.) interaction in the Bose-Hubbard model ()2.14p corresponds to A = +oo, 
or Amin = +00 for the IRH Bose-Hubbard model (j2.15p . This formally discards from the boson 
Fock space 5^s(A) all vectors with more than one particle at one site. 

Let $0 denote the vacuum vector in ^^^(A). Then the subspace ^^^•'^■(A) C 5^b(A), which 
corresponds to the hard-core restrictions, is spanned by the orthonormal vectors 

$^ = J] a: $0 , X C A . (3.1) 

Since the subspace 5^^'^ (A) is closed, there is orthogonal projection Pa onto ^^^■^■(A) such that 

tB'iA) = PdB{A) , (3.2) 
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and we get the representation 



(3.3) 



where the orthogonal comphment (^^^^'^■(A))-'- := (/ — P)^b{A). 

Since our main Theorem 12.11 is vahd for any Amin > and the estimate ()2.46|) is uniform in 
X'^, we can extend this theorem to the hard-core case by taking the hmit Amin +oo. 

For simphcity we consider the case of a sequence of non-random identical and increasing 
positive {X'^ = Xs > 0}^^ such that A^ — » +oo. 

Lemma 3.1 Let Xg +oo. Then for all ( ^ C : Im{Q ^ 0, and for any u ^ Q and u E C 
we have the strong resolvent convergence of Hamiltonians \2.2'J^) : 



hm {Hl{s,v)-CiyH = P 



, * G dB{A) , 
(3.4) 



where 



HX{s, v) := Ta + XsJ2 - 1) + J] e>x - Wiyao + 



van 



The same is true for approximating Hamiltonians \2.2J^ : 
hm {Er'"'\s,z,v)-CirH = 



p 



(3.5) 



(3.6) 



-1 



for any z E C and \I' G ^b{A). Here 



H 



u,appr ^ 



V\zf-\^(zao + zal)+NA + X,J2 



^x \J^X ' 



a;eA 



(3.7) 



Proof: By estimate ()2.32j) and ()3.5j) for < A^ < A^+i we get: 



V ^ 



A« 



uV) N^-V < H%{s, u) < m(s + 1, z/) . 



(3.^ 



So, for any uj E Vt and G C Hamiltonians ()3.5|) form an increasing sequence of self-adjoint 
operators, semi-bounded from below. Let {/i^ (z^, A)[\l/] := (^, Hj({s, z/)^f)5g(A)}^i be the corre- 
sponding monotonic sequence of closed symmetric quadratic forms with domains dom /i^(z^. A). 
Put 

g := Pi dom /i^(z/,A) , (3.9) 

and let Qq = Q he the closure of Q in the Hilbert space 5^b(A). Since for any u E Q and z^ G C 

hm {^,HX{s,u)<i/h,^A) = +oc , ^EitB^iA))^ , (3.10) 

As— >+oo 

one gets Qo = 5^^'^ (A) and the strong resolvent convergence ()3.4|) of Hamiltonians, see e.g. jDj, 
Ch.4.4 or |NZj . Lemma 2.10. (Note that for hard cores the space ^^'^'{A) is finite-dimensional, 
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which makes these arguments even simpler.) The strong resolvent convergence ()3.4|) of Hamil- 
tonians implies also 

hm ($,ff;;:(s,i.)$)5,(A)= (3.11) 

(<I>,P[Ta + 5^£X- v^(z7ao + z/a*)]^'^)5-= (A) , $ e ^?|-^-(A) . 

The same line of reasoning leads to ()3.6|) for approximating Hamiltonians. □ 
By the Trotter approximating theorem |RSlj the convergence fj3.4|) and ()3.6|) yields the 
strong convergence of the Gibbs semigroups: 

Corollary 3.1 The following strong limits exist: 

s- hm e-^«X(s,-) = e-^«i^.c.,AM , (3.12) 

where 

HL.,Ai^) ■■= P[Ta + E^>- - ^(^^0 + ^^o)]^ , (3.13) 

and similarly 

s - hm e-^^A"''"^^'^'^) = e-^^-T(^'^) , dom R^^ ^(u) = 5^-^-(A) , (3.14) 

As^+oo ' 

where 

K;cT(^' ^) •= 1^1' - ^(^«o + zal) + + 1)^- - ^(^^0 + iyal)]P , (3.15) 

with dom H';^^]^{z, V) = ds'-W- 

Since |e^'^*-^A'(*''^)~/^^A) j^^^ jg sequence of trace-class operators from Ci(5^b(A)) monotonously 
decreasing to the trace-class operator 

e-/3(^^jr....AM-/^A^A) eCi(^?^-^-(A)) , 

the convergence (j3.12j) can be lifted to the trace- norm topology, see [Zj. The same is true for 
fj3.14p . It then follows that the pressures also converge: 

Lemma 3.2 

hm p[HXis,u)]=p[H^,,,^^iu)] , (3.16) 
hm p[H-^{s, u)] = p[K:^J{z, u)] . (3.17) 

As— >+oo ' 

Since the estimate (j2.46j) is uniform in A > Amin > 0, we can take the limit A^ — > +oo to obtain 

< P [HL.Ai'^)] - p [H-::j{zA,M, u)] < 0{1/V) , (3.18) 

for all G fi, any /3 > 0, fi & M} and < Cq. Then, by the same line of reasoning as after 
fj2.46p in Theorem 12. H we obtain the thermodynamic limit of the pressure for the hard-core 
bosons: 
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Corollary 3.2 The pressure of the Infinite-Range-Hopping hard-core Bose-Huhhard model with 
randomness is given by 

PhAP,f^)= (3.19) 
sup {-r' + riE{lnTr(;j....)^ exp(/5P [(/i - - l)n, + r{a: + a,.)] P)}] , 

r>0 >- V « ' J 

cf. expression \2.19j} for finite A . 

Remark 3.1 To calculate the Tr over note that the boson creation and annihilation op- 
erators are quite different from operators : c* := Pa*^P, '■= Pa^P restricted to dom c* = 
dom Cx = ^b'^'; which occur in m^.iy\} . The major difference consists in their commutation 
relations: 



*\2 



[c.,c;] = 0, (xy^y), (c.)^ = (c:.) 
Taking the XY representation of relations \'J.2IJ^) : 







* I * J 



(3.20) 





1 



r 



1 




\3.19^) gives explicitly 



sup 

r>0 



E -(/i-e--l)+/3-Mn 



2cosh(l^v/(^737^^Tj^ 



(3.21) 



the grand- canonical pressure for the random IRH hard-core Bose-Hubbard model. 



3.2 Limit of Perfect Bosons 

The limit A — > is more delicate. For simplicity, below we assume that emin = 0. Then Hamil- 
tonian ()2.15|) for perfect bosons A!j^ = is non-negative, i.e. the corresponding pressure exists 
in a finite volume only for negative chemical potentials. There is an analogue of Lemma f3. 11 if 
we subtract from this Hamiltonian a term with fi < and assume u small enough: 

Lemma 3.3 Assume that e^i^ = and let A., \ 0. Then for yU < 0, for all C, E C : Im[Q ^ 0, 
and for any uj G VL, we have the strong resolvent convergence of Hamiltonians \2.2'J\) : 

\im{HX{s,u) - fiNA-Ciy^"^ = (3.22) 

for u E C, if |z^p < The same is true for approximating Hamiltonians \2.2J^ : 

lim {Hl''^\s, z, u) - fiN^ - C/)-'^ = (3.23) 
{V \zf - Wizao + zd*o) + J2i^x + 1 - f^)nx - Vv(udo + ual) - C/}"'^ , 

for any z eC, C eC: Im{C) ^ and ^ e dsi^). 
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Proof: The bound ()2.32|) now yields: 

HXis, z/, /i) := HXis, u) - fiNA > (-/i - \u\^)Na - V , (3.24) 

so that for |z/p + /x < 0, the operators {H'^{s, u, Ai)}^>i are positive. As in Lemma ITTl for these 
operators we define the corresponding closed symmetric quadratic forms by {h'^{u, fi, A)[^] : = 
(\1/, if^(s, z/, /i)\l/)5g(A)}^i. Note that they are monotonously decreasing and bounded from 
below, which implies that for any u E fl, u E C and A the operators {if^(s, i^,fi)}g>i converge 
in the strong resolvent sense, see e.g. |Kj, Ch.VIII, to a positive self-adjoint operator _ff^Q(z/, /i). 
Let us define the symmetric form 

h^J^]= lim/z-[$], (3.25) 

s— >oo 

with domain 

dom(0 = U domm. 

S>1 

It is known, [K] Ch.VIII, that if the form ()3.25j) is closable, then operator H'^q{u, fi) is associated 
with the closure h"^. By explicit expression of /i^(i/, /i. A) one gets that the limit form ()3.25p is 
closable (and even closed), since it is associated with the self-adjoint operator H'^{s = oo, z/, fi). 
Then the operator H^q{v,^) associated with the closure of ()3.25j) simply coincides with 
Ht{s = oo,i/,/i): 

h^J^ = ($ , [Ta + J](£- - ^)n,. - WiVao + K)] ^) , 

that proves ()3.22j) . 

A similar argument applies for the approximating Hamiltonians ()2.24j) . But, in contrast 
to the case of sources |z/p < that we can choose as small as we want to apply the main 
Theorem 12.11 the value of z will be defined by variational principle ()2.19p with A!^ > 0. Now 
the semi-boundedness of {H'^'"'^^^ {s, z, from below follows from the estimate 

I 1 2 

J2i^x + 1 - l^hx - WiiV + z)do + {u + z)al) >-V ^^^^ . (3.26) 

The rest of the arguments is identical to those for the operators ()3.24j) . or equivalently for the 
sequence {Hj({s, z^)}^>i , and goes through verbatim to give the proof of the limit ()3.23j) with 

H';^'''PP'{s = oo, z, u) := H';^';'""~{z, u). □ 

Corollary 3.3 In a full analogy with Corollary \S. 1\ and Lemma lS.^ the Trotter approximation 
theorem and the monotonicity of the operator families {H'^{s,i')}^y-^ , {H'^'"'^^^ [s , z , u)} ^^-^ 
yield 

\\m p[HXis,u)]=p[Hl,iu)] , (3.27) 

As — *U 

\im p[H-^'^^^^{s, z, u)] = p[H-:riz, u)] . (3.28) 

Notice that, similarly to the Weakly Imperfect Bose-Gas |ZBj . the estimate i2.4(^ for fi < is 
still uniform in X > 0. Therefore, we can take there the limit — > to obtain 

< P [HXM - P Ko''^(^A,4^), ^)] < 0{1/V) , (3.29) 
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for all u; G any P > and \u\ < —ji. Then, following the same line of reasoning as after 
\2.4t^ in Theorem \2. il we obtain the thermodynamic limit of the pressure for the perfect bosons: 



Po{P,fi< 0) = 

sup {-r^ + /3-^E{lnTr(5^)^ exp(/3 [(// - £^ - l)n, + rial + a,)])}} 



r>0 



cf. expression h2.1f^) for finite \, where all values of fi are allowed. Since we put e-^ 
variational principle in iS. .VCJj) implies: 

Po(/5, /i < 0) = r'E{lnTr(5,)^ exp{/3 [(^ - £^ - l)nj)} = 
/5-iE{ln[l-exp{/5(/x-e--l)}]-^} . 



(3.30) 

0, the 
(3.31) 



The convexity of {p [Hj( q{i' = 0)] }^ anc? i/ie thermodynamic limit po{P , fi) as the functions of 
< 0, together with the Griffith lemma, see e.g. |ZBj . yield the convergence of derivative with 
respect of n, i.e. the formula for the total particle density: 



p{P,H< 0) = E 



o/3(l+e"-At) _ I 



(3.32) 



Remark 3.2 As usually in the case of the perfect boson gas one recovers the value of thermo- 
dynamic parameters at extreme point fi = by continuation: fi —0: 



PoiP, /i = 0) := r {In [1 - exp{/3(-£^ - l)}]"^} 

1 



p(A/i = 0) :=E 



g/3(l+e-) 



(3.33) 
(3.34) 



In particular by \3. 34\ ) it gets clear that the gap (= 1) in the one-particle spectrum of the perfect 
boson gas Ta and emin = imply that the critical density 



PciP) := sup p{p, /i) = p(/5, /i = 0) 

/i<0 



(3.35) 



is finite, cf. i2.12^} and H2.1!^} . This opens a room for the zero-mode Bose condensation in the 
case of the random potential {s'^}^. 



4 Phase Diagram 

Here we analyse only the case, when e'^ is random, but the interaction couplings A|Jf = A > 
are fixed. 

To proceed we recall first the formulae determining the critical temperature /3c(p, A)~^ for 
the nonrandom case = 0. To this end we define, cf ()2.5()j) . 

p{p,p,\;r) := ^ In Tr^ exp(-/3 [/i„(/i. A) - r{a* + a)]) , (4.1) 

where 

hn{p, X) ■■= - P)n + Xn{n - 1) . (4.2) 
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Due to |BDj it is known that the critical temperature (and the critical chemical potential 
/ic(p, A)) are defined, as functions of the total particle density p, by two equations: 



Here 



and 



/(/3,/i,A;0) = 2, p 



/(/5,/i,A;0) 



^ oo 



ne 



-f3hn{tJ.,\) 



Zo(/5,/i, A)^ hn-iin, X) - hnin, X) 



(4.3) 



(4.4) 



Zo{f3,fi,X) = J2e~^'" 



n=0 



If 7^ and A > 0, then by the main Theorem O (see (EUHl), (^3^1 and (|01)) to obtain 
the equations for the critical temperature and the critical chemical potential we have to replace 
/i in ()4.3p by /i — ei^ and to average over e"^. This gives, instead of ()4.3p . the (gap) equation: 



E[p"iP,p-e^,X;0)] = 2 



and equation for density: 



p = E 



E 



ne 



-/3h„{/i-£",A) 



Zo(/5,/i-£-,A) ^ 
The case of A = is more subtle, and we begin with it the next subsection. 



(4.5) 



(4.6) 



4.1 Perfect bosons: A = 

Without loss of generality, we can assume that the random e'^ takes values in the interval [0, e]. 
In that case the maximal allowed value for n (i.e. the critical value) is still /i^ = 0, and the 
critical inverse temperature Pc '■= Pc{p, X = 0) is given (see ()3.34|) . ()3.35|) ) by: 



P 



E 



1 



g/3,{l+e-) 



(4.7) 



Remark that, irrespective of the distribution of e^, the equation ()4.7|) implies that the 

resulting Pc is lower than In ^1 + ^j, which corresponds to the nonrandom case e'^ = 0, i.e. 

disorder enhances Bose-Einstein condensation. We shall see (Sect. 4. 3. 3) that this is no longer 
true when A > 0, and even that the opposite holds, if A is small enough! 

Notice that formula ()4.7j) is in agreement with the general expression found in |L-Zj : 



P 



dM{E) 

pl3cE _ 1 



(4.8) 



where M{E) is the integrated density of states given by 



Ar(E) = a.s.- hm -#{z: Et<E}. 

v^oo V 



(4.9) 
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Here {E!^}-^-^ are the eigenvalues of the one-particle Hamiltonian with a random potential 

:= {tAu){x) + ^5Xx) , X e A, ue f)(A), (4.10) 

for i.r. kinetic-energy hopping, see 1^2.1^ . ()2.1()|1 . and ^{i : i^^^ < E} counting the number of 
the corresponding eigenfunctions (including the multiplicity) . It is known that for any ergodic 
random potential {£x}xeA^ limit ()4.9|) exists almost surely (a.s.) and that it is non-random, 
see e.g.jPF]. A contact between formulae ()4.7|) and ()4.8|) gives the following 

Lemma 4.1 The integrated density of states is equal to 

J\f{E) = P [e"^ < ^ - 1] = E [e{E - (1 + e"^))] . (4.11) 

Proof: For simplicity we consider the case of a Bernoulli random potential such that 

e'^ = e with probability p and e!^ = with probability 1 — p. (The proof of the general case 
is similar, but slightly more complicated.) In this special case, the right-hand side of ()4.1ip 
equals 

( 1 if E > 1 + e, 
F[e'^<E-l] = l 1-p iil < E <l + e, (4.12) 
[ if E < 1. 

Clearly, all eigenvalues {E'^}i>i of the Hamiltonian ()4.10p belong to the interval [1, 1 + e]. 
Since dim(f)(A)) = V, one gets ^{E) = 1, if E > 1 + e. Similarly, Af{E) = 0, if E < 1. 

Now suppose that £'€[1,1 + e). Since {£x}x£A Bernoulli random field, for given 

6 > 0, there exists c > such that with probability Pr > 1 — 6 the number of sites x G A with 
e'^ = e is in the interval {pV — c\/V,pV + cVV). Given a configuration for which this is the 
case, let A^ C A be the set where e't^. = e. Consider the states (p G f)(A) such that (f){x) = 0, if 
X ^ Ag and J2xeA '^(^) ^ ^- Then 

ihl<P){x) = - J2{<P{x) - <P{y)) + £^0(x) = ie + l)0(x) , x e A,. 

y=l 

The space of such eigenfunctions has dimension lA^I — 1, so that 

#{Er >E}> (|A,i - 1). 

Since (#{£;^ < E}) + {i^{E^ > E}) = V", for ^ oo we get 

^f{E) <l-p. 

Similarly, considering the eigenfunctions with supports concentrated on A^ = A \ A^ we obtain 

^f{E) >l-p. 

Together with ()4.12|) these estimates give the proof of ()4.11|) . □ 
The relations ()4.11|) show that the formulae ()4.7p and ()4.8|) are equivalent. For details of a 
general statement see e.g. jPFj Ch.II.5 . 
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4.2 Discrete random potential and A > 

We now consider the case with interaction A > 0, and first assume that the probabihty distri- 
bution of £^ is discrete. 

A particularly simple case corresponds to the hard-core boson limit A = +00, see Section 
3. Then by ()3.21|) the equations for the critical value of the inverse temperature Pc '■= Pc{p) = 
/3c{p, X = +00) for a given density p, reduce to the system: 



E 



tanhp{fi-e^ - l)/2 



1 



and 



P 



1 
2 



-E 



tanh -e'^ -1] 



(4.13) 



(4.14) 



The last equation ()4.14|) implies that for the hard-core interaction the total particle density has 
the estimate: p < 1. 



4.2.1 Bernoulli random potential in the hard-core limit A 



-00. 



A special case of a discrete distribution is the Bernoulli distribution, where e'^ = e with proba- 
bility p and e'^ = with probability 1 — p. We first consider the case A = +00. The equations 
(imSI) and (HTT4I1 then read, 



Fp,e{P = PcP) ■=P 



tanh ^Pdp 



fi — e — 1 



p) 



tanh \(3c{p — 1) 
fi — 1 



and 



1 

2 



ptanh -Pc{p ~ ^ ~ ly 



p)tanh-(3c{p - 1) 



P- 



(4.15) 



(4.16) 



Here a new phenomenon occurs for density p = 1 — p. To see this, we consider first a 
particular case of p = 1/2. Then p = 1/2, and by ()4.16|) we obtain, that the only possible 
solution for the corresponding chemical potential is p{p = 1/2) := p{p = 1/2, A = +00) = 
1 + 5/2. Inserting this value of p into (|4.15|) we get for the critical temperature: 

tanh = —e . 
4 2 

This equation obviously has no solution for e > 2. Therefore, there is no Bose-Einstein con- 
densation for Bernoulli random potential. Hp = p = 1/2, and e is greater than some critical 
value: Ecr = 2. 

One can check that the same phenomenon occurs for p 1/2 and for densities p = 1 — p, if 
e is large enough, but now the reasoning is more delicate. First of all, by (j4.15|) and tanhu < u 
we see that in any case there is a lower bound on the inverse critical temperature: 

Pc > 2. (4.17) 
Now assume that p < 1/2, i.e. p > 1/2. From ()4.16|) it then follows that for any e one has 

1 



0< p 



(4.18) 
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Indeed, if we suppose that < /U — 1 < e/2, then tanh |/5c(/U — 1) < tanh ^/5c(l + e — fi) and 
hence, by fl4.16|) . we get 

2p-l = ptanh^/3e(/i-£- 1) + (1 -p)tanh^/5c(/i- 1) 
< (1 - 2p) tanh ^pc{e + 1 - /i) < 1 - 2p , 

contradicting our assumption p = 1 — p, if /3c exists and is finite. 
Now notice that (|4.16p with p = 1 — p is equivalent to 



(4.19) 



1 — tanh ^Pc{£ + 1 — A^) 1 — p 
1 — tanh l(3c{f^ — 1) p 

The left-hand side of (I4.19p can be estimated from below as 

1 - tanh l/3,(e + 1 - /i) _ e/^-(/-'^-/^) + e"^--/^ M.-i-sm 
1 - tanh - 1) e-Mf^~^~e/2) ^ e 

Together with ()4.17|) this yield an upper bound for ()4.18|) : 

1 ll-pll-pl-2p 
< II -1 - -e < —In ^<-ln -< -. 4.20 

^ 2 pc P - 2 p 2p ^ ' 

But fl4.2()|) implies that ()4.15j) has no solution (3c, since for large e we obtain 

tanh - £ - 1) tanh - 1) 

P ^ + 1-P < 4.21 

p — e — 1 p — 1 

p 1 — p p 1 — p 

< —, r-^- r-, 1 r-^ < 1 • 



e + l-p p-1 e/2-{l-2p)/2p s/2 
We assumed that p < 1/2. Therefore by (j4.2H) . our conclusion is true, in fact, for 

e>l/p>2 = ecr. (4.22) 

The same result follows in the case p > 1/2, if we interchange p and 1—p and p — 1 and 1 + e — p 
in the above argument. 

Next we show that for any other p G (0, 1), i.e. for any p ^ 1—p, the critical Pdp) < +oo, 
i.e. for these densities one always has the Bose-Einstein condensation at low temperatures. 

To this end suppose that there is p* G (0, 1) such that p* 1 — p, but limp^p* Pdp) = +oo. 
Then the left-hand side of fj4.15p converges to 

lim Fp,,(/3, p) = Mp{p, e) := ^— + . (4.23) 

/3^oo \p — e — 1\ \p — l\ 

The number of solutions of equation ()4.15|) in the limit lim^^p. Pc{p) = +oo depends on the 
value of 5 > 0, but two singular points p = 1 and p = 1 + e of the function ()4.23|) ensure 
(for nontrivial values of the probability: p ^ and p ^ 1) that there are always at least two 
solutions: pi{e) < 1 and p2{e) > 1 + e of equation 

M^{p,e) = l. (4.24) 
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If linip^p. Pc{p) = +00, then for these two cases the equation ()4.16|) imphes: 

p* = hm Gp^siPcip), fJ'iie) = , 

p* = hm Gp^eiPcip),P2{£) = 1 • 
p^p* 

This contradicts our assumptions on p* and makes impossible the hypothesis hm^^p. Pcip) = 
+00. 

Notice that the function Mp{p,e) has a minimum 'Jl{e) G (1,1 + e). If Mp{jl{e),e) < 1 
(which is equivalent to e > Sp := l + 2^Jp{l — p)), then equation ()4.24j) has two complementary 
solutions /i=p(£:): 

P^{e) = '-^-PT ^(^)'-p(l^ , (4.25) 

such that 

1 < p,^{e) < Jl{e) < p+{e) < 1+ e . 
If limp^p. Pc{p) = +00, then for these two solutions equation ()4.1fij) implies: 

p* = Mm Gp,e{l3c{p),P^{e)) = 1 - p , 

This again contradicts our assumption about p*, and thus proves the assertion: Pdp) < +00 
for any p ^ 1 — p. 

Notice that by (|4.25p the equation Mp(/l(e),£) = 1 has a unique solution e = Ep < Ecr = 2, 
and one obtains Mp(jl{e),e) > 1 for all e < Ep, which excludes complementary solutions Pzf{e). 
On the other hand, if 

e > Ecr = max Ep = £p=i/2 , (4.26) 
p 

there are always complementary solutions ()4.25p . This may restrict the values of p, for which 
we have bounded critical Pdp), to a certain domain of densities. 

To this end we consider first the p-independent equation ()4.15p . Notice that Fp^^{P,p) is a 
monotonously increasing function of (3, so there is a unique solution Pdp) of equation (j4.15|) 
for a given p, if there is one. 

Since {ta.nh.u)/u < 1, then the left-hand side of ()4.15p is less than 1, for /3 < 2. On the 
other hand, as /5 — > 00, the left-hand side of ()4.15p converges to Mp{p,E). Since the function 
()4.23p is singular at p = 1 and p = 1 + E, a. solution 2 < Pdp) < +00 for a certain p always 
exists, and the set of those p is defined by the condition: 

Sp,, := {peR': lim Fp,,(/5, p) = Mp{p, e) > 1} (4.27) 

By (|4.23|) the set (j4.27p for e > is a compact in R^. If there are no complementary solutions 
p=p(£:), this compact is connected, but if 

E> Ecr . (4.28) 

it contains two domains separated by a gap: 

I{e,p) := (/i-(e), p+ie)), 
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see fl4.25|) . The gap I{e,p) C (1, 1 + e). There is no solutions /3c(/i) for /i e I{e,p) and for 



/i < (e + l)/2 - ^{(e-l)/2f-e{l-p) , 

or for 

/i > (5 + 3)/2 +^{{8 + l)/2f - 8(1 - p) . 

Hence, for large 8 ()4.28j) the set Sp^e is a union of two (separated by the gap I{8,p)) bounded 
domains, which are vicinities of singular points /i = 1 and fi = 1 + 8 . is in fact not the 

To understand, how the gap in the chemical potential for solution /3c(/x) modify the be- 
haviour of Pc{p), we have to consider the p -dependent equation ()4.16p . Notice that from ()4.16p 
one obtains Pc{p,p) as a function of two variables. Therefore, Pdp) is a solution of equation: 

/3,(/i) = /3,(/i,p) , (4.29) 

which in fact connects fi and p: JI^p), i.e. Pdp) = Pcipip)) = Pc{'p{p)yP)- 

Clearly, the left-hand side Gp^siP, p) is increasing in p and it tends to as — > — oo and to 
1 as /i — > +00. Excluding p = or 1, there is therefore a unique solution p{P,p) of ()4.16|) for 
each value of /?. As /3 — > 0, Gp^siP^p) tends to 1/2 at constant p. Therefore, if p 7^ 1/2 

lim u(3, p) = ±00 , 

13^0 

depending on whether p>l/2orp<l/2. 

On the other hand, in the limit /3 — > 00, we have that Gp^^iPjp)'- (c^) tends to 0, if p < 1; 
(6) to (1 - p)/2, if p = 1; (c) to 1 - p, if 1 < p < 1 + 8; (d) to 1 - p/2, if p = 1 + e, and (e) to 
1, if p > 1 + e. 

The (a) — (e) give relation between p and p for large (3 : if < p < 1 — p, we must have 
p(/3, p) — > 1 and, if 1 — p < p < 1, we obtain p(/3, p) — 1 + e, for /3 ^ 00. At p = 1 — p, we 
have to use the representation ()4.19|1 . that yields 

p(/5, p = 1 - p) = 1 + - ^ In ^ + oip-') , (4.30) 

if f3 is large. In particular, this justifies the remark ()4.22|) above about = 2, since 1 + e/2 
lies in the gap I{8,p) only if e > 2 = see (|4.25|) . 

Hence, it follows that for p 7^ 1 — p two functions of p corresponding to solutions (|4.29|) 
of equations (|4.15p . (j4.16p must intersect. On the other hand, (|4.22|) proves that they can not 
intersect for p = l— p,if£:> In fact, we can derive upper bounds for /3c(p) in the case 
p ^ 1 — p and |p — 1 + p| small. 

To this end we first consider the case p > 1 — p. Let us assume p < 1/2. (The case p> 1/2 
can be studied similarly.) Writing p=l— p + 6/2we present the equation ()4.1(ij) in the form 

ptanh^/3c(e + 1 -p) = (1 - p) tanh ^/3c(p - 1) + 2p - 1 - 6 . (4.31) 
Identity 1)4.311) implies that p > 1 +8/2, since otherwise we get a contradiction: 
l-2p + 6 = -ptanh ^/5c(5 + 1 - p) + (1 - p) tanh ^(3c{p - 1) < 
— ptanh -Pci^ + 1 — A^) + (1 — p) tanh -Pd^ + 1 — p) < 1 — 2p . 
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On the other hand, for e > 1, one gets the upper hmit fi < e + 1. Indeed, if we suppose the 
opposite: fi > e + 1, then (|4.16|) and the general fact that /3c > 2 (see (j4.17|) ) yield 

1 — 2p + 6 = p tanh -/3c(/i — £ — 1) + (1 — p) tanh -/3c(/i — 1) 
> {I — p) tanh ^/3c(/^ — 1) > (1 — p) tanhe. 
But this is impossible for (large) e verifying: 

Therefore, we obtain for /i the lower and upper bounds: 

1 + £/2 < /i < 1 +e . (4.33) 

Now identity ()4.3H) . together with the bounds ()4.33|) . inequality tanh(M) > 1 — 2e^^" and 
/3c > 2 (see ()4.17|) ). yields the estimates: 

A 9 1 A 

1 e^" < tanh -^Ae + 1 - /i) < 1 - -. (4.34) 

p p 2 p 

and hence, 

f3c< ^^-^ln(2p/5). (4.35) 

p — — 2e ^ 

The upper bound ()4.35p holds for example, if 5 < p/2 and e > ln(4/p). 

Now we consider the case p < 1 — p and suppose p < 1/2, since p > 1/2 can be studied 
similarly. Then we write: p = 1 — p — 6/2. Equation ()4.16|) now reads as 

(1 - p) tanh ^/3c(/i - 1) = ptanh ^/3c(l +e-^) + l-2p-6. (4.36) 

An argument similar to the case p > 1 — p shows that 

1 < p<l+e , (4.37) 
if e is large enough and S < 1 — p. Indeed, if we suppose the opposite: fi > I + e, then 

1 -2p-(5 > (1 - p) tanh ^/?c(/i- 1) > (1 -p)tanh£: , 
which is impossible for 

1, 2-3p-5 
e > - In — . 

2 p + S 

Similarly, if we suppose that fi < I, then ()4.36|) implies 

> ptanh ^/3c(l + e-p) + l-2p-6> ptanhe + {1 - 2p - 6) , 
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which is impossible if 5 < 1 — 2p , or if 1 — 2p < 5 < 1 — p and 

1 3p-l + 5 

2 1 — p — d 

Now, (jOF^ and (jOTjl imply that 



1 5 
tanh - 1) < 1 - —— . (4.38) 

In the case /U > 1 + |e this yields immediately the upper bound : 

/3.<?lnH(i^. (4.39) 

e 

On the other hand, ifl</i<l + £:/2, then by ()4.36|) and /9c > 2 we obtain 

1 1 
(1 - p) tanh -/3c(/i - 1) > ptanh -/^c^ + 1 - 2p - 5 > 

ptanh + 1 - 2p- S > p{l - 2e"^) + l-2p-5 = l- p- (5 - 2pe"^ . (4.40) 
Taking into account equation ()4.15|) and estimates ()4.38p . ()4.40|) . we get 

1 > > (3 c 



ln(2(l-p)/(5) ' 
that gives the upper bound: 

1 — p — — 2pe d 
4.2.2 Bernoulli random potential for the case A < +oo. 

We assume in this subsection that A > 5 + 1. If the repulsion is very large (A ^ e + 1), the 
analysis for p < 1 is then almost the same as above for A = +oo, whereas for p > 1, which is 
possible only for finite A, one needs some more arguments. 

Here we start with the estimate the first-order correction in A^^ to the value of £:cr(A = 
+oo) = 2. With this accuracy the equations (j4.5|) and (j4.6p can be approximated correspond- 
ingly by 

/ tanhl/3(p-.-l) 1 e-Z^g^-^)/^ A 

V 2A + £ + l-p coshi/3(l + e-/i)y ^' ^ 



/tanhi/3(/i- 1) 1 e^^'^-i)/^ 



V yu-l 2A + 1 -/i coshi/5(/i- 1^ 

and by ()4.16p as above. 

To see this, note that if p < 1, the dominant contribution in ()4.6p must come from the 
n = 1 term, i.e. we must have /ii < /;.2, so p < 1 + 2A + e. The other terms in (j4.6|) are then 
exponentially small and can be neglected, which leads again to (j4.16|) . 
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Now, because of the presence of e ^^'^ in the n = 2 term of ()4.4|) . it cannot be neglected in 
()4.5|) and we obtain: 

2» r e~^(l+^-^') - 1 g-/3(l+e-At) 

^ ' ■ 2- 



1 + e-/3(i+^-/^) [ /i-l-e l + 2A + £- /i 
2(1 -p) re-^(^-'^)-l e-^(i-^) 



1 + e-/3(i-M) [ ^ _ 1 1 + 2A - /i 



which is the same as (I4.42p . 

Similar to (j4.23|) the gap equation for 1 < fi < 1 + e can be obtained from (j4.42|) in the 
limit f3 oo: 

^ + (1 - (-^ + WT-^:^] = 1- (4.43) 



e + l — fi \fi—l 2A + 1 — /i 

If p = 1 — p, then by (j4.16p and (j4.3U|) we again obtain the limit: fi —>■ 1 + for f] ^ oo 
Inserting this limit into (j4.43|) we obtain 

2 , 2(1 -p) 



2A-ie 



1 . (4.44) 



Hence, by the reasoning similar to those after ()4.30|) . we obtain the critical value of the Bernoulli 
random potential ecr(A) the expression: 

^-^^) ^ l.^l_p)/X = 2 + 2(1 - P)/A + . . . , (4.45) 

which takes into account that A is large but finite. 

Another observation, which is related to the finiteness of A, concerns the value /5c(p = 1)- 
For hard-core bosons the arguments in the Sect. 4. 2.1 show that this value is infinite and the 
corresponding values of the chemical potential must be greater than 1 + e, see (I4.6|l . Now for 
finite A and fi > 1 + e the limit of ()4.42|1 . when (3 ^ oo, reads as: 

P ( ^ + TTT-^, ) + (1 -P) + WT-^ ] = 1- (4-46) 



2\ + l+e- fij 2\ + l-fi 

If p > 1, then we need to reconsider the density equation ()4.6|) . which has the form: 

p = v ^"=^ + {l-v) ^"=^ (4 47) 

Notice that ii (3 ^ +oo, then by ()4.2j) and ()4.47p one obtains the following limits: p — > 1, 
when p G (1 + e, 1 + 2A) , p — 2 — p, when p G (1 + 2A, 1 + 2A + e), and p — 2, when 
p G (1 + 2A + £, 1 +4A). 

Therefore, at p = 1 for large (3 we can ignore in ()4.47|) the terms higher than see fl4.2|) . 
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and write in this limit: 



P 



1 + g-/3(l+£-At) _^ g-2/3(l+A+e-^j) 



1 + 2e"'^(^+ 2A+£-m) 



1 + g-/3(At-l-e) _^ g-/3(l+2A+e-M) 

1 + 2e-^(i+2^-^) 



1 + g-/3(M-l) + g-/3(l+2A-M) 



(4.48) 



1 + p (e 



'/3(l+2A+e-M) 



-/3(At-l-£) 



) 



v) 



-/3(1+2A-m) 



)■ 



This yields 



,2/3/. 



32/3(l+A) 



1 — p + pe 

1 — p + pe 



0e 



-Pe 



_P_g2/3(l+A+i£)_ 



1 — p 



e as 



The chemical potential defined by equation (j4.47p therefore tends (for p = 1) to 1 + A 

Therefore, inserting this into ()4.46p we obtain the estimate for the value of repulsion Ac,i 
that ensures that /?c(p = 1) = +oo in the presence of the random Bernoulli potential: 



KAe) = - 



9 + 2e{l - 2p + -e) 



(4.49) 



Remark 4.1 In the absence of disorder, i.e. if e = 0, the critical value of lambda is Ac,i = 3 
as opposed to Ai = |(3 + v^) as suggested in \BDf . The reason is the same as above for Scr, 
namely, the graph of fi{j3,p) at p = 1 tends to 1 + X as j3 +oo and this lies in the gap only 
if X>3. Similarly, the next critical values are given by 



Kkie = 0) = 2k + 1. 



(4.50) 



Remark 4.2 In Sect. 4. 2.1 we notice a new phenomenon specific for the random case: diver- 
gence of (3c at p = 1 —p for hard-core bosons, cf. Figure 1 for p = 1/2. Instead of fixing X, 
fixing e > 2 it follows from \4.44 ) that there is a critical value of the repulsion Ac,i-p(e) (instead 
of e as in (|^.^5| )j so that Pdp = 1 — p) diverges for X > Xc^i^p{e) in the presence of the random 
Bernoulli potential: 

Ki~pie) = I + ^T^l^- (4.51) 
This critical value is not evident from Figure 1 as e = 2. 



Remark 4.3 In Sect. 4. 1 we remarked that the critical temperature for free bosons increases 
due to disorder. We also remarked that for the interacting case this is a more subtle matter, 
since it depends on the value of repulsion. For large repulsions close to e.g. Ac,i(£ = 0) = 3, 
we get by that 

p^(p = 1; A = 3, £ > 0) < Pc{p = 1; A = 3, e = 0) = +00 . (4.52) 
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This lowering of j3c{p = 1) can be explained intuitively as follows. At density p = I, there is 
one particle per site. If e = there is a penalty for a particle to jump to an already occupied 
site, so the preferred state is where the particles are at fixed sites, which is almost an eigenstate 
of the number operators Ux for each site. This prevents Bose condensation. (This argument 
was presented also in IBDf . ) However, if e > 0, then the lattice splits into two parts with 
energies and e, and a particle jumping from a site with energy e to a site with energy loses 
an energy e, which counteracts the gain of A. This creates more freedom of movement and 
therefore promotes Bose condensation. On the other hand, for a fractional value of the p in the 
neighbourhood of p = I — p, the critical temperature decreases with increasing e as can be seen 
from Figure 1. 

Now consider the case p > 1. From equation ()4.47|) we see that at fixed p G (1,2 — p), 
/i — > 1 + 2A and for p G (2 — p, 2), — > 1 + 2A + e as /3 — > oo. 

For the case p = 2 — p, we have to expand ()4.47p . as above for p = 1, see ()4.48|) . but to take 
into account that /i G (1 + 2A, 1 + 2A + e): 

f 1 + 2e"'^(^+ 2\+e-^l) ^ 



1 + g-/3(M-l-£) _^ g-/3(l+2A+£-M) 

g/3(l+2A-M) I 2 



1 + g/3(l+2A-M) + g-2/3(M-l-A) 
~ 2-p + p (^g-/3(l+£+2A-/i) _ g-/3(M-l~e)^ _ 
(1 _ p)g-/3(/..-l-2A) _ 2(1 _ p)g-2/3(M-l-A)_ 

This yields that e-^^^^-^-^A) ^ g_/3(i+e+2A-M)p/(i _ p) for large /5, i.e. /i ^ 1 + 2A + if 
p = 2 — p and /3 — > oo. 

For /i !^ 1 + 2A + l^:, one has hi{^ — A) < /i2(p — A). So that the p-terms in ()4.42p 
are unchanged, but A) > /i2(yU, A) < h^i^fi, A), if A > e/4, which corresponds to our initial 
hypothesis about the value of repulsion: A > 1 + e. Hence, the (1— p)-terms are now dominated 
for large (3 hj n = 2 and (j4.42j) read as 

P f g-/3(l+£-M) _ I g-/3(l+£-M) 



1 + g-/3(i+e-/.) 1^ fi-l-e l + 2\ + e~fi 

1 _ p r g-2/3(l-M+A) _ g-/3(l-/.) g-2/3(l-M+A) ^ 

+ g-/3(l-M) + g-2/3(l-M+A) Y /i - 1 -2A ^ ^ 1 + 4A - /X J ^' 

In the limit /3 — >^ oo we obtain from this relation the gap equation 

p ( i + ) + (4.54) 

^'-^^ C-L2A + 1 + 4LJ='- 
Inserting p. = 1 + 2X + into (|4.54p leads to 

-e^ - (2A - 1 + 2p)e + 8A = 0. (4.55) 

Solutions of (|4.55|) are: 

^ (2A - 1 + 2p) ± v/(2A - 1 + 2p)2 - 16A . (4.56) 
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Hence, there is a solution that for large A has the form: 

4?nA) =4(1 + 1^)+... , (4.57) 

or other way around, for a given e we have: 

K,=2-,{e) = ^|f^ . (4.58) 

Clearly, this critical value only applies if e > 4 and p > 1/2. The top graph of Figure 1 
illustrates this behaviour at p = 1.5 for e = 4.5 and A = 10. 

The critical /3c{p) for the Bernoulli distribution with p = 1/2 and 5 = 2 is shown in Figure 
1 for a number of values of A. Notice in particular that e < ecr(A), see (|4.45|) . for all finite A, 
so that i3c{p = 1/2) < +00. 

Also, for A = 3.3, one obtains Pdp = 1) < +00 because 3.3 < Ac,i(£ = 2) = (3 + vT3)/2, 
see dOni). 



6 
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2 
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0.5 1 1.5 2 

Figure 1: Pc as a function of the density p in the case of aver- 
aging over two energies: and e = 2 with equal probabilities, 
for various values of X: A = 3,3.3,4,6,10 and +00. The top 
graph corresponds to the case e = 4.5 and A = 10. 




4.2.3 Trinomial distribution: A = +00. 

We also briefly consider the trinomial distribution, taking for simplicity equal probabilities, i.e. 

Pr = 1/3 

Pr = 1/3 (4.59) 
Pr = 1/3 . 

For hard-core bosons, A = +00, equation ()4.13|1 for the critical value of Pdp) takes the form: 



tanh |/3(/i — 1) tanh |/3(/i — 1 — le) tanh |/3(/i — 1 — e 



/i 
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+ 



p — 1 — e 



1 . 



(4.60) 



The density equation ()4.14|) now reads as 



p = - + -( tanh - 1) + tanh 



2 ' 



tanh 



(4.61) 



Then by the same analysis as in Sect. 4. 2.1 one gets from ()4.6H) : 



hm p(/3, /i) 



if p < 1 
1/6 if fi = l 

1/3 if 1 < I2<1+e/2 

1/2 if /i = l+£/2 

2/3 if 1 + e/2 < /i < 1 

5/6 if /i = l+ 5 

1 if /i > 1 + e . 



Other way around this can be also expressed as: 

' 1 

l + e/4 
l + £/2 
1 + 3e/4 
l+e 



lim 



if < p < 1/3; 

if p=l/3; 

if 1/3 < p < 2/3; 

if p = 2/3; 

if p > 2/3. 



Again, similar to the reasoning in Sect. 4. 2.1, the inserting ofp=l + £:/4or/i=l + 3£:/4 into 
the limiting equation ()4.60|) for j3 — > +oo yields the critical value of the random potential: 

^cr = y . (4.62) 

Therefore, (similar to the Bernoulli case for p = 1/2) the condensation of hard-core bosons is 
absent at densities p = 1/3 and p = 2/3, if £ > This phenomenon of course persists for 
A < +00 and there are similar suppressions of Bose condensation at p = 4/3,5/3, etc., if e is 
large enough. 



4.2.4 Trinomial distribution: A < +oo. 

For A < +00 there is a similar enhancement of Bose condensation at p = 1 as for the Bernoulli 
distribution, but the effect is stronger. This can be seen in Figure 2. The explanation is similar 
to that in Remark l4.3^ except now the lattice splits into 3 equal parts with energies 0, e/2 and 
e. Particles can jump from a singly-occupied site with energy e to a singly-occupied site with 
energy or e/2, thus compensating for the energy penalty of A due to double occupation. 

By equation ()4.14|) for ()4.59j) we obtain that at p = 1, p(/5, p) — * 1 + A + e/2 as — > +oo. 
The gap equation ()4.60|) then reduces to 

1 1 1 _ 

A - e/2 ^ A ^ A + e/2 ~ ' 

We can solve it for e provided A > 3: 

e^X) = 2X^^ . (4.63) 

Thus, Bose condensation is absent, if A > 3 and e < ecr(A). 

Figure 2 shows /9c(p) for a fixed e = 10 and for values of A > 3. Then e > ecr(A = 3, 4, 6), 
but e < ecr(A = 8) = 13.52, which excludes condensation at p = 1 in the latter case. 
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Figure 2: (3c as a function of the density p in the case of a 
trinomial distribution with width e = 10 for A = 3, 4, 6 and 8. 



4.2.5 General discrete distribution. 

The same phenomena persist for higher numbers of random potential energy values, but the 
critical value ecr(A) becomes rapidly very large. Figure 3 shows the case of a distribution with 
equals probabilities Pr = 1/10 at 10 equidistant values of e'^ (with maximal value e = 10) 
for A = 8. Clearly, condensation is suppressed at p = 1/10, . . . ,9/10 and p = 1,2 but not at 
corresponding fractional values above 1, cf. Figure 2. 
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Figure 3: (3c as a function of the density p in the case of 
averaging over 10 energy values with width e = 10 for A = 8. 
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4.3 Continuous distribution 



4.3.1 The case A = +00. 

Consider a random potential with homogeneous distribution between and e. In case A = +00 
the equations ()4.13|) and ()4.14|1 become 

1 rtanhi/5(/.-l-x)^^^^ 



e Jq /i — i — X 
and 

1 r 1 

-/ tanh-/3{i2-l-x)dx = 2p-l . (4.65) 
^ Jo 2 

The latter has sense only for < p < 1 and can be solved exactly for /j,: 

^ g/3(/.-l-£)/2 + g-/3(^-l-£)/2 = 2P - 1 , 

and hence 

1 1 sinh ^/3ioe 
,(ftp) = l + 5.+ ^ln-5^^^^. (4.66) 

As /3 ^ +00, the expression ()4.66|) takes the form 

lim /i(/?, p) := 7J(p) = l+ep, 0<p<l, (4.67) 

,3^+00 

whereas p(p = 0) G (—00, 1] and 7l(p = 1) G [1 + e, +00) for extreme values of density, i.e., the 
inverse function is 

( p < 1 

p(p) = < (p-l)/e l<p<l+e (4.68) 
[1 1 + e <fi . 

Then by ()4.64p and ()4.68p we obtain for p = 1 in the limit (3 +00: 

1 = - dx , 

e Jq p - 1 - X 

or we get explicitly the value of the chemical potential 

p(p=l) = l + -^^>l+£, 

i — e 

and similarly 

p p = = 1 - — - < 1 • 

1 — e ^ 

Hence, for hard-core bosons the critical /3c(p) is infinite at extreme densities p = 0, 1 for any 
value e > of the uniform continuous distribution. 

If < p < 1, then solution of the equation ()4.65p in the limit P — > +00 is ()4.67p . whereas the 
integral in ()4.64|) diverges. Therefore, if the critical /3c(0 < p < 1) exist, it must be bounded. 
Moreover, since (ta.nh.u)/u < 1, by ()4.64|) we get for it a bound from below: 2 < /3c(0 < p < 1). 

To prove the existence and uniqueness of /5c(0 < p < 1) consider first ()4.65|) for p < |. Then 
by virtue of ()4.66j) for any finite /? the solution fi{P,p) increases from —00 to 1 + e/2 when p 
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changes from to 1/2. For this variation of chemical potential the integral in the left-hand side 
of (j4.64p increases monotonously from to its maximal value given by 

1 r tanhi/3(a;-e/2) 
fi=l + e/2) = - / '"^^ ' U x . (4.69) 

^ Jo X — Ej I 

Indeed, 

1 /tanhi/5(/i-l) tanhi/?(^-l-£)\ 

9^/ /?, /i = - — \ — > 

/i— 1 /i — 1 — £ / 

for /i < 1 + 5/2. The integral in (|4.69p is obviously an increasing function of /3. So, there exist 
/?o > 2 such that the maximal value of integral /(/5o,/i = 1+^/2) > 1. Hence, for any (3 > (5o 
there is a unique density < p(/?) < 1/2 such that 

/(/5,/i(/5,p(/5)) = l . (4.70) 

Notice that by ()4.66|) p) is increasing of the both arguments: (3 and < p < 1/2. Hence, 
to satisfy ()4.70|1 p(/3) must be decreasing function of /3, i.e., the inverse function (3^ = Pc{p) is 
also a decreasing with limp^o Pc{p) = +00 and limp^i/2 Pdp) > Po- 

Similar arguments are valid for 1/2 < p < 1. Whereas p,{f3,p) is still increasing function of 
p, the integral I{(3, p) now decreases with p from its maximal value ()4.69p to 0. Therefore, I3c = 
Pc{p) is a monotonously increasing function of p with limp„,i/2 /3c(p) > Po and limp^i /?c(p) = 
+00 , i.e. with a minimum at p = 1/2 as we have seen for discrete distributions and hard-core 
bosons. 

4.3.2 The case of large A < +00. 

By virtue of equations (j4.5|l and (j4.6|l . for A < -|-oo, the Bose condensate is still suppressed at 
p = k. 

The analysis is very similar to the case e = 0. In the limit /3 — >■ -|-oo by ()4.6p the density 
tends to {k = 0,1, . . .) 

if p < 1 

p(p, (3)^ { A; + i(p - 1 - 2A;A) if 1 + 2A;A < /i < 1 + 2A;A + e 

k + 1 if l + 2A;A + £<p< 1 + 2(A; + 1)A. 

(To see this note that if 1 -|- 2kX < p < 1 + 2kX + e then the term e~^^''+^ dominates for 
X < p — 1 — 2kX and the term e"^'^* dominates for x > p — 1 — 2kX.) Clearly, if < p < 1 then 
for solution of (j4.6|) one gets as above: p{[3, p) 1 + pe when [3 +00. If p = 1, we need to 
approximate ()4.6p more carefully: 

1 e/3(M-l-a;) _|_ 2e2/3(At-l-a;-A) 

-dx 



5 1 + ^|3{^l~l-x) g2^(/i-l-x-A) 



Working out the integral, we find that p(/?, p=l)^l + X + }^eas[3^ +00. More generally, 
ii p = k, p{[3, p = k) 1 + {2k — 1)A -|- ^e. For large (3, the gap equation (j4.5|) becomes 

k k+l . , 



p-l-2{k-l)X)-x l + 2kX + x + 2X-p 
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Inserting /i = 1 + {2k — 1)A + we obtain that 



k 



X + — X 

This gives for the critical values of repulsion: 



+ 



A 



dx = 1. 



1 e^/(^^+^) + 1 



(4.71) 



It is easy to see that this is larger than for non-random case Ac,fc(0) = 2k + 1 and agrees with 
the value mentioned above at e = 0, see Sect. 4. 2. 2 . 

Figure 3 shows the phase diagram for A = 10 with e = 3, taking an average over a uniform 
distribution corresponding to 10 equidistant random values of e'^ in the interval [0, 3]. It shows 
that this already approximates the continuous case quite well. 
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Figure 4: j3c as a function of the density p in the case of a 
near- continuous distribution: averaging over 10 energy values 
with width e = 3 for A = 10. The lower graph is the case 
without randomness. 



4.3.3 The case of small A > 0. 

We finally consider the case of small A. Figure 4 shows that, in contradistinction to the case 
A = 0, for small A the critical (3c{\,s) > I3c{\ = 0,e = 0), i.e. it is larger than that at e = 0! 

This can be understood as follows. Whereas in the free case A = 0, we must have < 0, 
when A > 0, this is no longer so. In the limit A — ^ 0, we can replace e~^'^"'^^'^'^ in the expression 
(|4.4p for p"{P, jj, A; 0) occurring in the gap equation (j4.5p by e^^^~^\ Replacing also — /i„ 
(see (|4.2j) ) by /i — 1 the series (j4.4j) can be summed and we obtain for (j4.5|) : 

- / dx = 1 . 

e Jq 1 + X - h 

If e = this leads to the free gas critical value fi = 0, but for £ > we obtain 

e" -1 - e 

,= '-^^>0. (4.72) 
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Figure 5: (3^ as a function of the density p in the case of 
averaging over two energies and width e = 2 for small A = 
0.1. For comparison, the lower graph shows the case without 
randomness. 



Similarly, the density equation ()4.6|) now reads as 




g/3(l-/^+a;) 



dx = p . 



(4.73) 



By ()4.72|) we can approximate for small e fi hj fi ^ e/2 and inserting it in ()4.73p we find 



V P/ 

Notice that this argument also applies in the case of a discrete distribution, see Figure 5. 

5 Conclusion 

We conclude by few remarks concerning our results and open problems. Summarizing the most 
striking observations about the model considered in this paper, we have seen that at large values 
of the on-site repulsion with a discrete distribution of the random single-site particle potential, 
the disorder causes a suppression of Bose-Einstein condensation at fractional values of the 
density. On the other hand, the suppression of Bose-Einstein condensation at integer values of 
the density observed in the absence of disorder is lifted. For continuous distributions we found 
that the critical temperature decreases with increasing disorder for non-integer densities. 

We have have concentrated here on the case of uniformly distributed random external po- 
tential. Nonuniform distributions as well as a random on-site interaction may also be of interest 




(4.74) 



By convexity of the function (e'^'^^"'"^) 



1) \ we conclude for solution of the equation ()4.74|) that 
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and give rise to new phenomena. Of course, all our results concern the infinite-range-hopping 
model. It would be of considerable interest to extend our results to the short-range hopping 
model. 
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